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Further Observations of the Tilted Planet XO-3: A New Determination 
of Spin-Orbit Misalignment, and Limits on Differential Rotation* 
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Abstract 

We report on observations of the Rossiter-McLaughHn (RM) effect for the XO-3 exoplanetary system. 
The RM effect for the system was previously measured by two different groups, but their results were 
statistically inconsistent. To obtain a decisive result we observed two full transits of X0-3b with the 
Subaru 8.2-m telescope. By modeling these data with a new and more accurate analytic formula for the 
RM effect, we find the projected spin-orbit angle to be A = 37.3° ±3.0°, in good agreement with the 
previous finding by Winn et al. (2009). In addition, an offset of ^ 22 m s~^ was observed between the 
two transit datasets. This offset could be a signal of a third body in the XO-3 system, a possibility that 
should be checked with future observations. We also attempt to search for a possible signature of the 
stellar differential rotation in the RM data for the first time, and put weak upper limits on the differential 
rotation parameters. 

Key words: stars: planetary systems: individual (XO-3) — stars: rotation — techniques: radial 
velocities — techniques: spectroscopic — 



1. Introduction 

Since the discovery of the first transiting planet, many 
groups have been studying the stellar obliquities (spin- 
orbit angles) of planet-hosting stars through measure- 
ments of the Rossiter-McLaughlin (RM) effect. The RM 
effect is a distortion of stellar spectral lines that occurs 
during transits, originating from the partial occultation of 
the rotating stellar surface. It is often manifested as a pat- 
tern of anomalous radial velocities (RVs) during a plane- 
tary transit (Queloz et al. 2000; Ohta et al. 2005; Winn 
et al. 2005; Narita et al. 2007; Triaud et al. 2010). By 
modeling the RM effect, one can determine the angle A 
between the sky projections of the stellar rotational axis 
and the orbital axis. The statistics of the spin-orbit angle 
A should provide a clue to the formation and evolution 
of close-in giant planets (hot-Jupiters and hot-Neptunes). 
Since 2008, many transiting systems with significant spin- 
orbit misalignments have been reported (e.g. Hcbrard 
et al. 2008; Narita et al. 2009b; Pont et al. 2010). This 
has attracted much attention to the importance of dy- 
namical mechanisms for producing close-in planets, as 
well as tidal evolution of planets and their host stars 
(Fabrycky & Tremaine 2007; Wu et al. 2007; Nagasawa et 
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al. 2008; Chatterjee et al. 2008; Triaud et al. 2010; Winn 
et al. 2010). 

In this paper, we present the measurement of the RM 
effect for the XO-3 system. The XO-3 system was discov- 
ered by Johns-KruU et al. (2008). Photometric follow-ups 
by Winn et al. (2008) allowed the system parameters to 
be refined. The large mass of the planet {Mp = 11.79 ± 
0.59Afjup) and its eccentric orbit (e = 0.260 ±0.017) at- 
tracted further interest in this system. 

Hcbrard et al. (2008) detected the RM effect with the 
SOPHIE instrument on the 1.93 m telescope at Haute 
Provence Observatory (OHP). They found A = 70° ± 15°, 
suggesting a significant spin-orbit misalignment for the 
first time among the known planetary systems. On 
the other hand, Winn et al. (2009) independently mea- 
sured the RM effect with the High Resolution Echelle 
Spectrometer (HIRES) installed on the Keck I telescope, 
and found A = 37.3° ±3.7°, which differs by more than 2 a 
from the former result. 

The reason for the discrepancy was unclear, but it may 
indicate the presence of unknown systematic errors in one 
of the datasets, or even in both. It is equally possible that 
the discrepancy should be ascribed to the different tech- 
niques adopted in modeling the RM effect. Hcbrard et al. 
(2008) and Winn et al. (2009) both used analytic formu- 
lae to compute the anomalous RVs, but the former group 
used a formula that was based on a calculation of the first 
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moment of the distorted line profile, while the latter group 
used a formula that was calibrated by numerical analysis 
of simulated RM spectra. Recently Hirano et al. (2011) 
presented a new and more accurate analytic formula for 
the RM effect, showing in particular that the RM velocity 
anomaly depends on many factors such as the rotational 
velocity of the star, the macroturbulent velocity, and even 
the instrumental profile (IP) of the spectrograph, not all 
of which were considered in the previous literatures. 

Specifically, for rapidly rotating stars like XO-3, the 
velocity anomaly calculated by Hirano et al. (2011) dif- 
fers strongly from the simpler, previous analytic descrip- 
tions based on the first-moment approach (Ohta et al. 
2005; Gimenez 2006). When the incorrect relation is used 
between the RM velocity anomaly and the position of the 
planet, the results for A may be biased. 

In order to resolve the disagreement, and obtain a de- 
cisive result for the angle A with fewer systematic errors, 
we observed another two full transits of X0-3b with the 
High Dispersion Spectrograph (HDS) on the Subaru 8.2- 
m telescope. We also applied the new analytic formula by 
Hirano et al. (2011) to model the RM effect with greater 
accuracy. We find that the best-fit value for A based on 
our new measurements is very close to that reported by 
Winn et al. (2009). 

We describe the detail of the observation in Section 2. 
The data analysis procedure and the derived parameters 
are presented in Section 3. Section 4 discusses the compar- 
ison with the previous results, and considers the possible 
effect of the stellar differential rotation. 

2. Observations 

We observed two complete transits of X0-3b with 
Subaru/HDS on November 29, 2009 and February 4, 2010 
(UT). We also obtained several out-of-transit spectra on 
each of those two nights as well as on January 15, 2010 
(UT). The out-of-transit spectra were obtained in order 
to help establish the Keplerian orbital parameters of the 
system. We adopted a typical exposure time as 600-750 
seconds, and chose the slit width as 0.4", corresponding 
to the spectral resolution of ^ 90,000. We used the Iodine 
cell for precise RV calibration. 

We reduced the images to one-dimensional (ID) spectra 
using standard IRAF procedures. The typical signal-to- 
noise ratio was ^100 per pixel in the ID spectra. We 
then processed the reduced spectra with the RV analysis 
routines for Subaru/HDS developed by Sato et al. (2002). 
Table 1 gives the resulting RVs (corrected for the motion 
of the Earth) and the associated errors, which are com- 
puted from the dispersion of RVs that were determined 
from individual 4 A segments of the spectrum (Sato et al. 
2002). We obtained a typical RV precision of 11-14 m s^^. 

3. Analysis and Results 

3.1. Fit to the Subaru RV data alone 

We determined the projected spin-orbit angle A in sev- 
eral steps. First, in order to provide an independent de- 
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Table 1. Radial velocities measured with Subaru/HDS. 

Time [BJD (TDB)] Relative RV [m s~^] Error [m s~^ 



termination, we use only the transit data from our new 
Subaru observations. Since those data alone are insuffi- 
cient to determine all the Keplerian orbital parameters of 
the system, we also use the out-of-transit RV data points 
from OHP/SOPHIE (Hebrard et al. 2008). This essen- 
tially provides an independent determination of A since we 
do not use the in-transit RV data from OHP/SOPHIE. 

Our model for the RVs is similar in some respects 
to the previous analyses by Narita et al. (2009a) and 
Narita et al. (2010). Each RV data set (Subaru/HDS and 
OHP/SOPHIE) is modeled as 

Kiodoi = K[cos{f + w) + ecos(n7)] + Aurm + Toffsct, (1) 



Z^OOIU^. ( UOl ( ^ 




1 4 


0A^^^ (\A 71 1 R1 zL 




±0.0 


0A^^^ (\A 71 Q^({A 




1 4 


Z'^iOoiD'y:. ( z ( o\j^ 




1 1 

10. 1 


OA^^^f{A 7'X^C\^A 

Z'^zOOiD'y:. ( OOUO^^: 


404 1 


11.0 


Z^rJrJXU^. ( 4:ZOU4: 


ouo.y 


lO.O 


9zL'^^1 (\A 7^C\^AA 


ooo.z 


19 8 
iz .o 


9zL'^^l (\A 7^R97^ 

Z^OOIU^. ( JOZ ( O 


977 8 


1 
lO.U 




9'^9 


19 7 

IZ. ( 


0A^^^ (\A 77*^74^ 

Z'4:001U'4:. t ( O i ^0 


1 QS 
lyo.o 


1 9 ^1 
IZ . o 


0A^^^ fi/L 7R1 AR^ 


1 7 


1 1 

10.1 


9zL'^^l f{A 78090^ 


1 40 4 


1 9 Q 
iz . y 


9zL'^^l f{A 7Qf{QA^ 


1111 
i i J L 


19 8 
iz .o 


9zL'^^l f{A RCiAf\7^ 


119 9 


19 9 
iz . z 


9zL'^^1 f\A Rl 940^ 


1 no 7 




lO.U 




1 "^n Q 


1 1 Q 
11. y 


OA^^^(^A R97R7^ 


1 ^7 7 


19 7 
IZ. i 


9A^^^ f{A k'\^f{f\^ 

Z'4:001U'4:.00<JUUO 


1 94 n 


19 4 

IZ .4: 


9A^^^ f{A k'\Q^^^ 


1 fifi 4 


91 8 
Z 1 .o 


9zL'^'^91 1 79004^ 

Z^rJrJZXX. ( ZUU4:U 




14 ^1 
l^.O 


9zL'^'^91 1 790^1 

Z^rJrJZ X X . ( Zc/U lU 


774 n 


1 7 

lO. ( 


9zL'^'^91 1 7^^807^ 
z^oozii. ( ooy ( J 


S44 1 


14 9 

1^. z 


0A^^0'^^ 700/1/10 
z^oozoi. ( uy^^u 


4Q9 S 


14 


9A^^'y'^^ 71R91Q 
z^oozo-L. ( loziy 


^94 8 


lO.O 


9A^^'y'^^ 79^QzLQ 
Z'4:OOZO-L. ( zoy4:y 


uUU.O 


1 Q 
lo.y 


Z'4:OOZOl. ( oouoo 


480 9 


lO.O 


OA^^O'^1 7A^ A^ r 


497 R 
^Z i .0 


10.0 




444 9 


1 1 

lO.l 




41 fi 1 


19^ 
IZ . o 






1 8 
lo.o 


9A^^'y'^^ 779"^!^*^ 


98'^ fi 
zoo.u 


1 4 Q 
1^. y 


9A^^'y'^^ 7R000R 


ouo. 1 


14 9 

1^. z 




990 1 

ZZU. J- 


1 
lO.U 


9/lc:c;0Q1 ICi^^lA 
z400zoi. ( yoo ( 4 


iOO.D 


19 9 

iz.z 


2455231.803313 


179.2 


13.4 


2455231.811063 


132.0 


12.8 


2455231.818802 


127.6 


12.8 


2455231.826532 


70.5 


12.3 


2455231.834271 


72.9 


12.0 


2455231.842001 


122.6 


13.2 


2455231.849740 


139.6 


13.1 


2455231.857469 


110.5 


13.2 



No. 



where K is the orbital RV semi-amphtude, / is the true 
anomaly, e is the orbital eccentricity, w is the angle be- 
tween the direction of the pericentcr and the line of sight, 
and finally 7 is a constant offset for the data from a given 
spectrograph. 

The RM velocity anomaly Awrm is modeled with 
Equation (16) of Hirano et al. (2011). In order to compute 
Awrm , we adopt the following values for the basic spectro- 
scopic parameters; the macroturbulence dispersion C, = 6.0 
km s~^, the Gaussian dispersion (including the instrumen- 
tal profile) (3 = 3.0 km s~^, and the Lorentzian dispersion 
7 = 1.0 km s^^. These values are taken from Gray (2005) 
and from the comparison with the numerical simulations 
by Hirano et al. (2011). Also, we assume the quadratic 
limb darkening law with ui = 0.32 and U2 = 0.36 following 
Claret (2004). 

We fit the two RV data sets (Subaru and OHP) by min- 
imizing 
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where V^^^ is the observed RV value labeled by i while 
Knodei corresponds to Equation (1). The uncertainty for 
each RV point is expressed by ct^*^. Since we do not 
have any new photometric observations of the transit, 
we fix the photometrically measured parameters to be 
RpjRs = 0.09057, a/i?, = 7.07, and io = 84.2° from the 
refined parameter set by Winn et al. (2008). The remain- 
ing parameters are K, e, zu, 7offsct (for each data set), 
the rotational velocity of the star u sin is, and the spin- 
orbit angle A. We allow all the parameters to vary freely 
to minimize x^, using the AMOEBA algorithm. We add 
the stellar jitter of cTjittor = 13.4 m in quadrature to 
the RV uncertainties in Table 1 so that the reduced in 
the global RV fitting becomes unity (after adding an ad- 
ditional parameter to allow for an offset between the two 
Subaru transits, as explained below). This jitter is ac- 
counted for in estimating the uncertainty for the system 
parameters in Table 2. 

By fitting the Subaru/HDS data along with the out-of- 
transit OHP/SOPHIE data, we find the spin-orbit angle 
to be A = 36.7° ± 3.0°. This is in agreement with the pre- 
vious finding by Winn et al. (2009), and in disagreement 
with the previous finding by Hcbrard et al. (2008). The 
reduced chi-squared is = 1.14. Interestingly, when we 
plot the residuals between the Subaru/HDS data and the 
best-fit model, we find a small negative trend as a function 
of time over the 67-day span of the observations. To show 
this, we plot our new out- of -transit RV data as a function 
of BJD in the upper panel of Figure 1, along with the 
best-fit curve (red). The residuals from the best- fit curve 
are shown at the bottom. This trend cannot be corrobo- 
rated or refuted by the previously published observations; 
the RV precision obtained by Johns-KruU et al. (2008) 
and Hebrard et al. (2008) was insufficient, and the precise 
RV measurements of Winn et al. (2009) did not cover a 
sufficiently long observation period. 
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Fig. 1. (Upper) New RV data outside of transits, obtained 
with Subaru/HDS. (Lower) The orbit of XO-3b based on the 
measurements with Subaru/HDS (blue), and the previously 
published RVs obtained with OHP/SOPHIE (black). For this 
figure, a linear RV trend (7) was fitted to the data and then 
subtracted. For each of the figures above, the best-fit model 
is shown as a red curve and the RV residuals from the best-fit 
model are plotted at the bottom. 
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Table 2. The best-fit parameter sets. 



100 



Parameter 



(A) Subaru (B) Subaru + Keck 



K [m 

vsinis 
AH 
7 [m s 



[km s 



-1 day-i] 



1499.5 ± 9.9 
u.zooy^Q QQ27 
347.4 ±1.4 
17.0 ±1.2 
37.3 ±3.0 
-0.322 ±0.088 
0.91 



1494.0 ± 9.5 
0.2883 ±0.0025 
346.li}i 
18.4 ±0^8 
37.4 ±2.2 
-0.320 ±0.088 
1.00 (fixed) 



This RV trend miglit indicate a possible additional body 
in the XO-3 system, but it is obviously premature to con- 
clude so only with the 3 epochs of data. Future obser- 
vations are needed. For the present purpose, in order to 
account for the offset in the overall RV between the dif- 
ferent transit epochs, we introduced an additional model 
parameter 7, representing a constant radial acceleration. 
Wc then refitted the Subaru RV data. The results from 
this fit are given in column (A) in Table 2. The uncer- 
tainty for each parameter is derived by the criteria that 
Ax^ becomes unity. The inclusion of the constant ac- 
celeration improves the reduced chi-squared significantly 
(x^ = 0.91 from 1.14 in the absence of 7) and the best-fit 
RV acceleration is 7 = —0.322 ±0.088 m s^^ day~^, indi- 
cating a 3.6ct detection. Since the two transit observations 
are separated by 67 days, the RV offset between the two 
transits is estimated as ~ 22 m s~^. The resultant RVs 
as a function of the orbital phase are shown in the lower 
panel of Figure 1. 

3.2. Joint fit to the Subaru and Keck data 

Now that we have seen that our new results by 
Subaru/HDS support the previous RM measurement by 
Keck/HIRES (Winn et al. 2009), we would like to try to 
combine the two independent measurements (Subaru and 
Keck) and carry out a joint analysis in order to derive the 
parameters with greater precision. 

We fit all of the transit data from Subaru/HDS 
and Keck/HIRES, and also the out-of-transit data from 
OHP/SOPHIE. We allow for a constant RV acceleration 
7, as before. We estimate the best-fit values for K, e, 
zu, vsinis, A, and 7 as in Section 3.1. The results are 
summarized in the column (B) of Table 2. Most of the 
values are very close to the best- fit values in case (A). 
The projected rotation rate of usinis — 18.4 ±0.8 km s~^ 
is in good agreement with the spectroscopically measured 
value (wsinjs = 18.54 ±0.17 km s~^, Johns-KruU et al. 
2008). The resulting phase- folded RV anomalies during 
transits are plotted in Figure 2, in which the Keplerian 
motion and the linear RV trend are subtracted from the 
data. The RV data taken by Subaru/HDS are indicated 
in blue for the first transit and purple for the second tran- 
sit, and those by Keck/HIRES are shown in green. The 
red solid curve is the best-fit curve based on the analytic 
formula of Hirano et al. (2011). 
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Fig. 2. RV data spanning the transit, after subtracting the 
orbital contributions to the velocity variation, and also a lin- 
ear function of time. The plotted data includes the new 
Subaru/HDS data (blue for the fist transit on UT 2009 Nov. 
29 and purple for the second transit on UT 2010 Feb. 4) and 
the previously published Keck/HIRES data taken on UT 2009 
Feb. 2 (green). The RV residuals are plotted at the bottom. 

4. Discussion and Summary 

We have investigated the RM effect for the XO-3 sys- 
tem, which was the first confirmed system with a signifi- 
cant spin-orbit misalignment (Hebrard et al. 2008). The 
new spectroscopic measurements including two full tran- 
sits taken by Subaru/HDS and the new analysis method 
using the analytic formula for the RM effect by Hirano et 
al. (2011), found the spin-orbit angle of A = 37.3° ±3.0°, 
supporting the result by Winn et al. (2009) based on the 
measurement with Keck/HIRES. The joint analysis of all 
the RV data sets covering three transits with Subaru/HDS 
and Keck/HIRES have shown that the projected stellar 
spin velocity estimated by the RM analysis well agrees 
with the spectroscopically measured value. 

Our analysis also detected an RV trend, or at least RV 
offsets, among the three epochs of the Subaru/HDS obser- 
vations. The cause of the extra RV variation is not clear. 
It is possibly an indication of a third body in the sys- 
tem: a stellar companion (binary), or an additional mas- 
sive planet. Nevertheless it should be noted that this star 
is known to have a high "RV jitter" of around 15 m s"-'^, 
and the precise physical causes and timescales of the jitter 
are not known. It is possible for starspots or other surface 
inhomogeneities being carried around by stellar rotation 
to produce a systematic offset in RV observations con- 
ducted on a single night. Since the rotational velocity 
of the star is large for a planet-hosting star, even a rela- 
tively small spot could cause an apparent RV anomaly in 
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a similar manner as the RM effect. For example, the RV 
acceleration of 22 m s^^ could be caused by a very dark 
spot whose size is only 0.002 of the total stellar disk area. 
The best way to investigate these possibilities is with ad- 
ditional measurements of the out-of-transit RV variation, 
with a precision better than 15 m s~^. 

As for the results for A, we would like to understand the 
reason for the discrepancy between the OHP/SOPHIE rc- 
suhs, and the Subaru/HDS + Keck/HIRES rcsuhs. To 
this end we try several additional tests. As we have 
pointed out, Hcbrard ct al. (2008) employed the analytic 
formula based on the first moment of the distorted line 
profiles to describe the RM effect (Ohta et al. 2005). For 
rapidly rotating stars, however, the RM velocity anomaly 
computed from the first moment significantly deviates 
from that based on the cross-correlation method (Hirano 
et al. 2010). Therefore, we reanalyze the OHP data using 
the new analytic formula by Hirano et al. (2011) to see if 
the original estimate for the spin-orbit angle A was biased. 
Instead of fixing the stellar spin velocity wsini., as done 
by Hebrard et al. (2008), we allow it to be a free param- 
eter, and fit all the OHP RV data (Hebrard et al. 2008). 
The resulting spin-orbit angle is A = 58.8° ±8.9° and the 
stellar spin velocity of vsinig = 15. 9 ±2. 6 km s~^. The cen- 
tral value for A approaches our new results (37.4° ± 2.2°). 
but they still disagree with each other with > 2cr. This 
shows that a biased model played only a minor role in the 
discrepancy. The major reason seems to have been sys- 
tematic effects in the OHP/SOPHIE dataset, perhaps due 
to the short-term or long-term instrumental systematics 
(instability) for fainter objects as reported by Husnoo et 
al. (2011). 

Incidentally, with only a small modification, the ana- 
lytic formula by Hirano et al. (2011) can also be used to 
calculate the RM velocity anomaly in the presence of dif- 
ferential rotation (DR) . The detection of DR would be of 
great interest for understanding the convective/rotational 
dynamics of the host star. Furthermore, it allows a possi- 
bility to break the degeneracy between the projected and 
the real three-dimensional spin-orbit misalignment angle 
by inferring the inclination angle of the stellar spin axis 
with respect to the line of sight (ig), an angle that is or- 
dinarily not measurable with the RM observations (if the 
star is a solid rotator). 

Since XO-3 has a comparably large wsinis, our new data 
may provide a good opportunity to search for the signa- 
ture of DR, or at least to put constraints on the degree of 
DR quantitatively. 

To model DR, wc introduce two major parameters: the 
stellar inclination is and the coefficient of DR, a. The 
stellar angular velocity $7 as a function of the latitude I 
on the stellar surface is written as = Qcqi^ ~ asin^ I), 
where f2cq is the angular velocity at the equator (Reiners 
2003b) . We step through a two-dimensional grid in a and 
cos is, and for each grid point we fit the RVs with the six 
parameters listed in Table 2. We compute the resulting 
at each point {a, cosig). Wc note here that the DR 
of our Sun is well described by a ~ 0.2. This also seems 
to be a typical value of other stars based on the spectral 




Ax =1.0 



Ax = 0.0 
>^ 



-0.05 





a 



0.05 0.1 0.15 0.2 



Fig. 3. Contour plot of Ax'^ in the space of ttie DR parame- 
ters a and is. The confidence region where Ax^ < 1.0 is sur- 
rounded by the red solid curve. We also show the confidence 
boundary of Ax^ = 2.30 by the black dashed line, which de- 
termines the Icr region in a two-dimensional parameter space. 

line analysis of (Reiners & Schmitt 2003a), although those 
authors also point out that some stars may have "anti- 
solar" like differential rotations in which a < 0. Thus, our 
grid extends from —0.2 < a < 0.2 and < cos is < 0.95. 
The case where cos is > 0.95 is very unlikely because the 
star would need to be rotating unrealistically rapidly to 
give the observed value of wsinis. 

Figure 3 shows contours of Ax^ = ~ Xmin ™ ^^'^ ("^i 
cosis) plain. The location of the best-fit model (defining 
the condition Ax^ = 0.0) is plotted with a black cross. 
This figure shows that with the current RV data, we are 
only able to provide fairly weak constraints on the param- 
eters. We are able to rule out the far upper left and right 
corners of this parameter space, corresponding to Solar- 
like DR viewed at low inclinations. Wc can rule out much 
stronger levels of DR (|a| ^ 0.5) regardless of orientation, 
but such strong levels of differential rotation are unlikely 
in any case. 

The non-detection of DR may be ascribed to the large 
stellar jitter of the host star (15 m s~^). This is often 
typical of relatively hot and rapidly rotating stars such 
as XO-3. The best cases for studying DR through the 
RM effect would be somewhat cooler stars that are still 
moderately rapid rotators (« 5-10 km s~^), for which a 
greater signal-to-noise ratio can be obtained. 
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